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In order to describe the exotic nuclear structure in unstable odd-yl or odd-odd nuclei, the de- 
formed relativistic Hartree Bogoliubov theory in continuum has been extended to incorporate the 
blocking effect due to the odd nucleon. For a microscopic and self-consistent description of pairing 
correlations, continuum, deformation, blocking effects, and the extended spatial density distribution 
in exotic nuclei, the deformed relativistic Hartree Bogoliubov equations are solved in a Woods-Saxon 
basis in which the radial wave functions have a proper asymptotic behavior at large r. The for- 
malism and numerical details are provided. The code is checked by comparing the results with 
those of spherical relativistic continuum Hartree Bogoliubov theory in the nucleus ^^O. The prolate 
deformed nucleus ^^C is studied by examining the neutron levels and density distributions. 

PACS numbers: 21.60.-n, 21.10.-k, 21.60.Jz 



(N 
> 

m 

rn 

O 



X 



As one of the most fascinating phenomena found in 
exotic nuclei, the nuclear halo contains many new and 
interesting features, such as a largely extended spatial 
density distribution, dilute and pure neutron matter, the 
Borromean property, and the coupling between bound 
states and the continuum Q. In the past, several rela- 
tivistic Hartree Bogoliubov theories with or without the 
Fock term have been developed for a self-consistent de- 
scription of spherical halo nuclei [@-@|. In order to in- 
vestigate halo phenomena in deformed nuclei in a mi- 
croscopic and self-consistent way, a deformed relativistic 
Hartree Bogoliubov (RHB) theory in continuum has been 
developed for even-even nuclei 0"0|. Within this model 
a decoupling in shape between the core and the halo has 
been predicted in some deformed nuclei close to the neu- 
tron drip line, e.g, in '^'^''^^Mg |, §. 

In order to describe the exotic nuclear structure in un- 
stable odd- A or odd-odd nuclei, the blocking effect of one 
or several nucleons has to be taken into account. In the 
present Letter, the deformed relativistic Hartree Bogoli- 
ubov theory in continuum [0-^ is extended to incorpo- 
rate the blocking effect due to an odd nucleon. In such a 
way, pairing correlations, continuum, deformation, block- 
ing effects, and extended spatial density distributions in 
exotic, odd-A or odd-odd nuclei can be taken into ac- 
count microscopically and self-consistently. 

In order to treat pairing correlations, the quasi-particle 
concept is adopted and the ground state of an even-even 
nucleus |$) is represented as a vacuum with respect to 
quasi-particles [hol. 



/3fc|$o)=0, foranfc = l,...iV, 
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where N is the dimension of the quasi-particle space and 
the quasi-particle operators /3j^,/3j, are obtained by the 
Bogoliubov transformation from the particle operators 
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The matrix W of the coefficients U and V is unitary 
which guarantees that the quasi-particle operators satisfy 
the same anti-communication rules as the particle opera- 
tors [O. Starting from the bare vacuum |0), the ground 
state of a system with an even number of particles |$o) 
can be constructed as. 



i*o)=n^fcio)' 



(3) 



where fc = 1, . . . A'^. For an odd system, in practice, the 
corresponding ground state can be constructed as one 
quasi-particle state I'I'i), 



\^i) = Pl\^^) = ii\{ih\Q), 



(4) 



where /3{ corresponds to the quasi-particle state with the 
lowest quasi-particle energy. In other words, the one 
quasi-particle state |$i) is the vacuum with respect to 
the set of quasi-particle operators (/S^, . . . , /3^) with 



P[^Pl P'2=P2, ■•■: P'n=Pn, 



(5) 
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and the exchange of the operators fj^ \-r i^i 
set of quasi-particle operators {(3[, . . . , (3'^, /3^^, . . . , /3^), 
which corresponds to the exchange of the columns 
{Ui^,Vi^) < — > iVil,U*^) in the matrix W. That is, the 



blocking effect in the odd system can be realized by ex- 
changing the creator /3j with the corresponding annihila- 
tor f3^ in the quasi-particle space. Accordingly, the block- 
ing effect in a multi-quasi-particle configurations can be 
treated. 

The covariant density functional theory has been ap- 
plied to describe successfully nuclear structure over the 
entire periodic table [[11 
goliubov equations 13 



[LSI . The relativistic Hartree Bo- 
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for the nucleons read, 



d^j 



-A* 



A 
-ho + X 






= Ek 



Vk 



,(6) 



where Ek is the quasiparticle energy, A is the chemical 
potential which guarantees the proper average particle 
number and ho is the Dirac Hamiltonian, 



hoir, r')^a.-p + V{r) + I3[M + S[r)). 
The scalar and vector potentials 



(7) 



(8) 



S{r) = gaa{r), 

V{r) = g^Lu^ir) + gpTsp"{r) + e^^A"{r), (9) 

depend on the scalar field a as well as the time-like com- 
ponents uj'^, p^, and A^ of the iso-scalar vector field w, the 
3-component of iso- vector vector field p and the photon 
field A. 

For a fully paired and axially symmetric deformed sys- 
tem with the time reversal symmetry, the projection of 
the total angular momentum on the symmetry axis n is 
conserved and each single particle state has a degeneracy 
of two. The RHB equation (R) can be reduced to half 
dimension M ^ N/2 and can be decomposed into de- 
generate blocks with quantum numbers +Q or —Q. The 
corresponding density and abnormal density matrix have 
dimension M and read 



n = V*V 



(10) 

(11) 



where V and U are the coefficients in the corresponding 
subspace. 

For an odd system with the /cb-th level blocked in the 
+Vl subspace, the time reversal symmetry is violated and 
there appear currents in the system. These currents show 
an axial symmetry, i.e., Vi remains a good quantum num- 
ber, but the quasi-particle energies are no longer degener- 
ate for the two subspaces, because the subspace with -1-17 
contains the odd particle and the corresponding subspace 
with —17 contains an empty level. Therefore, in princi- 
ple, we have to diagonalize the RHB equation twice, one 
for the subspace with -l-fi and the other for the subspace 
with — f2. Using the equal filling approximation which is 
usually made |l^, |l7|, we neglect currents and average 
in a statistical manner over the two configurations of a 
particle in the +Vt space and a particle in the — $7 space. 
The corresponding currents cancel each other and in this 
way we obtain in each step of the iteration fields with 



TABLE I. Ground state properties of ^'"'O from deformed 
RHB and RCHB calculations. _Rn, -Rp and Rt refer to neutron, 
proton and total root mean square radii in unit of "fm" . i?part , 
E„, £5°", E^, EZ°'\ Ep, Seoul, £pair and Et^t refer to energies 
of particle, a field, non-linear term of a field, uj field, non- 
linear term of lu field, p field. Coulomb field, paring energy 
and total binding energy in unit of "MeV" . 
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2.83 


2.83 
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2.57 
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FIG. 1. (Color online) Neutron single particle levels in the 
canonical basis and their occupation probability v for ^' O 
from the deformed RHB and RCHB calculations. 



the time reversal symmetry. In practice we average the 
density matrix p and symmetrize the abnormal density 
K in two subspaces and replace the two densities by 

P' = Pm.m + \ {UkMl - ^kjl) ' (12) 



"^MXA/ ^ 2 V^^b^k, 



^LUk^) 



(13) 



where Vk^ and Uk^, arc column vectors in the matri- 
ces V and U corresponding to the blocked level. Note 
that Eqs. ([l^ ) and ( pli| ) are the same as those given in 
Refs. ]16|, [l7|| if one considers the time reversal symmetry 
and that the dimension of the densities here is M = N/2. 
In the following applications we use the density func- 
tional PKl |lq] and a zero range density dependent pair- 
ing force [Q, P^ The deformed RHB equations are solved 
in a spherical Dirac Woods-Saxon basis p9| and the de- 
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FIG. 2. (Color online) Neutron density profiles of ^^C and ^'^C 
in deformed RHB calculations with the parameter set PKl. 



formed potentials and densities are expanded in terms 
of the Lcgendrc polynomials Px{cos6) (for details see 
Refs. 0^)- In order to check the accuracy of the de- 
formed code with the blocking effect we investigate the 
spherical nucleus ^^O and compare the results of the de- 
formed RHB code allowing only the spherical compo- 
nents of the fields, i.e., A = 0, with those obtained in 
the spherical relativistic continuum Hartree Bogoliubov 
(RCHB) theory Q. The ground state properties from 
the deformed RHB and the RCHB calculations are given 
in Table B. It is clearly seen that both calculations agree 
well with each other. For example, the difference between 
the total binding energies is about 0.02 MeV, which cor- 
responds to an accuracy of about 0.01%. 

In order to check the details of the deformed RHB cal- 
culation, we show in Fig. m the neutron single particle 
levels in the canonical basis for ^^O in comparison with 
the spherical RCHB results. The length of each level is 
proportional to the occupation probability. In the spheri- 
cal RCHB calculation, the blocked orbital is ld5/2- In the 
deformed RHB code, a spherical solution is enforced by 
allowing only A = component in the Legendre expan- 
sion mentioned earlier and therefore the three sublevels 
of the ld5/2 orbital with fl"" = 1/2+, 3/2+, and 5/2+ arc 
degenerate. The present results in Fig. H are obtained 
with the level J7'^ = 1/2+ blocked in the deformed RHB 
code, but blocking any of these three levels gives the same 
results. From Fig. n one finds a good agreement between 
the results from the deformed RHB and RCHB calcula- 
tions. 

After checking the code in the spherical case we now 
turn to an application of the deformed RHB theory for 
an exotic and deformed nucleus with an odd number of 
particles. The deformed halo candidate nucleus ^^C ||2^ 
is studied with the parameter set PKl ||lj]. In contrary to 
the spherical even-even nucleus ^*C, the deformed RHB 
theory predicts a deformation /3 = 0.25 for the odd-A 
nucleus ^^C. The neutron root mean square radii of ^'^C 
and ^^C are 2.56 fm and 2.79 fm respectively. 

The deformed RHB calculations lead for ^^C to a de- 



FIG. 3. (Color online) Neutron single particle levels in the 
canonical basis and their occupation probability v'^ for ^^C in 
the deformed RHB calculation. 



formed density distributions with an axial symmetry. In 
Fig. the neutron densities along the symmetry axis 
Pn{'rz,r±^ ~ 0) and perpendicular to the symmetry axis 
Pnifz = 0, r_L) are plotted as dashed and dashed-dotted 
lines. We also include the spherical density for ^^C as 
a reference by a solid line. It is interesting to see that, 
in the direction perpendicular to the symmetry axis, the 
neutron densities are almost the same for ^''C and ^^C 
at least when r^ < 6 fm. Along the symmetry axis the 
neutron density of ^^C extends much further than that 
of ^^C. This is partly due to that ^^C is prolate and the 
weakly bound 1/2+ level is occupied. 

The single particle levels of ^^C in the canonical basis 
are plotted in Fig. 0. Since it is a deformed nucleus, the 
lP3/2 orbit is split into two levels with Q,^ = l/2~ and 
^^ = 3/2^ respectively. There is one neutron occupying 
in Q.^ = 1/2+ level near the threshold. The occupation 
probability v^ = 0.5, indicating that it is averaged over 
the two configurations with ^^ = ±1/2+. Because of 
the deformation, this level is a mixture of the spherical 
orbits ld5/2 (62%) and 2si/2 (36%). The weakly-bound 
feature and the relatively large s-wave component of this 
level results in that the neutron density of ^^C extends 
further along the symmetry axis. 

In summary, the blocking effect due to an odd nucleon 
is incorporated in a deformed relativistic Hartree Bogoli- 
ubov (RHB) theory in continuum in order to describe 
odd-A or odd-odd exotic nuclei. The formalism is briefly 
presented and the numerical checks are carried out by 
comparing the results of the deformed RHB code for the 
spherical nucleus ^^O with those obtained from the spher- 
ical RCHB code. As a first application, the deformed 
nucleus ^^C is studied by examining the neutron levels, 
deformation as well as density distributions along and 
perpendicular to the symmetry axis. Along the symme- 
try axis the neutron density of "'^^C extends much further 
because ^^C is prolate and the weakly bound valence level 
with ri'^ ~ 1/2+ has a relatively large s-wave component. 
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